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ABSTRACT 

We present the DONUTS autoguiding algorithm, designed to fix stellar positions 
at the sub-pixel level for high-cadence time-series photometry, which is also capable 
of autoguiding on defocused stars. DONUTS was designed to calculate guide correc- 
tions from a series of science images and re-centre telescope pointing between each 
exposure. The algorithm has the unique ability of calculating guide corrections from 
under-sampled to heavily defocused point spread functions. We present the case for 
why such an algorithm is important for high precision photometry and give our results 
from off and on-sky testing. We discuss the limitations of DONUTS and the facilities 
where it soon will be deployed. 

Subject headings: Astronomical Techniques 



1. Introduction 



Many areas of modern astronomy require evermore precise time-series photometry. Examples of 
such areas include: determining the rotation rates of stars, measuring irregularities in asteroid light 
curves, characterising stellar flares, pulsations and variability and the characterisation of transiting 
exoplanets and eclipsing binaries. In order to achieve high-precision, potentially sub-millimagnitude 
(mmag), photometry of bright (photon noise limited) stars, systematic noise in the data must be 
minimised. As tracking errors and mechanical flexure cause stars to drift around on a CCD they 
sample different pixels over time. Additional systematic noise may be introduced by uncorrected 
pixel-to-pixel sensitivity variations, errors from flat fielding, vignetting and changes in seeing. This 
potentially limits the achievable accuracy of the photometry. In un-guided observations such drifts 
are inevitable. Even in autoguided observations mechanical flexure and differences in atmospheric 
refraction between the autoguiding and science cameras will cause stellar images to drift on the 
CCD. 

If the positions of stars can be held fixed in place on the detector throughout a time series of 
observations, the effects of most of the systematic errors previously mentioned may be reduced or 
eliminated. However, this is not always possible. Theoretically the best poin t spread function (PSF) 
for accurate photometry is a flat topped, top hat-like shape (JHowelll l200d ) . Orthogonal Transfer 
(OT) CCDs have the ab ility to shift their collected charge during an integration to create such a 
PSF. lHowell et al.l (J2003l ) investigate the potential use of OT CCD s for high precision photq i netry . 



A more commonly used te chnique is to defocus the telescope (e.g. ICharbonneau et al.ll200d . 



Southworth et al. 



2009 



M 



1201 



2006 



Barros et al. 



2011 



Defocusing allows the bulk of the pixels within 
the PSF to remain the same even in the presence of slight drifts, hence diluting any pixel-to-pixel 
sensitivity variations. It also limits the effects of changes in seeing as the amount of defocus applied 
typically dominates the size and shape of the PSF. Finally, defocusing also allows for higher signal- 



to-noise per unit time observations of bright objects by allowing longer integration times without 
saturating, simultaneously reducing overheads and increasing efficiency of the duty cycle. 

We present a new autoguiding algorithm, DONUTS, which has been designed to fix telescope 
pointing to < 0.2 pixels for high-cadence time-series photometry. DONUTS calculates guide cor- 
rections directly from the series of science images in real time, essentially eliminating the effects of 
mechanical flexure and differential refraction between traditional off-axis autoguiders and science 
imagers. The algorithm is also designed to calculate guide corrections for a range of PSFs, from 
under-sampled (FWHM < 2 pixels) to heavily defocused (FWHM ~ 40 pixels), so as to maintain 
all the benefits of defocused photometry that were previously described. Such an algorithm is an 
essential requirement for the new ground-based transiting exoplanet survey Next Generation Tran- 
sit Survejo(NGTS) which aims to push ground-based transit detection limits to the mmag regime 
to routinely discover Neptune and super Earth-sized exoplanets around bright (V< 13) stars. 

Transit photometry is crucial in determining the physical properties of exoplanets. Obser- 
vations of multiple transits allow the orbital period of the planet to be tightly constrained, and 
modelling the transit's shape allows the orbital inclination and the planet's size relative to its host 
star to be deduced. When combined with radial velocity measurements of the host star's reflex 
orbital motion the sin (i) ambiguity on the planet's mass is removed and so its bulk properties 
can be calculated and compared with, or used to constrain, planetary formation and evolutionary 
models. 

This paper is organised as follows, in Section [2] we describe the DONUTS algorithm used to 
measure translational shifts in a series of science images. In Section [3] we present the results of our 
off-sky testing of DONUTS and in Section [5] we summarise th e results from on-sky testing at th e 



Near Infra-red Transiting Exoplanet Survey (NITES) telescope (jMcCormac et ahllAfcA submitted! ). 
Finally, in Section [4.11 we discuss the potential use of DONUTS, comment on its limitations and 
summarise the work presented here. 



2. The algorithm and simulated data 

The typical procedure followed by the DONUTS algorithm in order to provide autoguiding 
during observations is as follows. After acquiring a target a series of science exposures is started. 
When the first science image has been read out it is analysed using the method outlined below and 
becomes the reference image to which subsequent images of this field will be compared to determine 
guide corrections. When the second science image is acquired it is analysed in a similar manner 
and pixel shifts in the X and Y directions are calculated. These pixel shifts are then converted to 
telescope coordinates (e.g. right ascension and declination), sent to the mount (possibly through 
a control loop, e.g. proportional integral derivative) and the process is repeated for each image 
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acquired. 

For simplicity, in the remainder of this section we wih assume that a 40 x 40 pixel image 
subsection is being used for guiding on a star that has been defocused such that the secondary 
obstruction is visible in the PSF (hence the name DONUTS). As shown later in Section[3l DONUTS 
is capable of calculating guide corrections from a full CCD array in a fraction of a second, therefore 
using the positional information from all the stars in the field, while maintaining sub-pixel accuracy. 
Hence, we envisage that the full frame mode of operation will be used in practice. 

After slewing to a target, a reference image of the field is acquired and analysed by DONUTS. 
The image is collapsed along the X and Y directions to create two 1-D projections: 



Image = I{x,y), (1) 

X,,i=(x,yi{x,y)], (2) 



y / 

X / 



Y,,i=\y,}^I{x,y)\. (3) 

The projections, X^ef and Y^ef; are saved in memory and used as references for subsequent 
images of this field. On acquiring new images this process is repeated and two comparison projec- 
tions are created, Xcomp and ycomp- To determine the linear shifts in X and Y the corresponding 
reference and comparison projections are cross correlated using Fast Fourier Transforms (FFTs). 
This method was adopted here as the cross correlation of two sets of two 1-D projections was found 
to be faster than cross correlating two 2-D images. One of the key advantages of autoguiding with 
DONUTS is that we avoid determining the centroid of our individual guide stars which is difficult 
for both under sampled and heavily defocused PSFs. 

In X, we take the FFT of both the reference and comparison projections: 



FFTx,^, = FFr(X,ef), (4) 

FFTx^on., = FFT{X,,^^). (5) 

The complex conjugate of the FFT of the comparison projection is then multiplied by the FFT of 
the reference projection: 

$x = FFT* (Xeomp) X FFTx^^f . (6) 

We then take the Inverse Fast Fourier Transform (IFFT) (px of the previous cross correlation <I>x 
and search the array for the location of the maximum (/>Xmax • 



</>x = /FFT($x). (7) 

The location, locx, of (/'Xmax ^^ the array of (px is directly related to the shift in pixels between 
the two projections analysed. Figure [U shows example guide corrections calculated for simulated X 
projections of a slightly defocused (such that the simulated PSF becomes double peaked) guide star, 
both in the case where the overall intensity of the comparison star flux has increased uniformly and 
where the intensity of the PSF is varying slightly from pixel-to-pixel in the projection. Assuming 
our reference point as the middle of the CCD array (i.e. pixel 20) for the reference star (see Fig [1] 
solid red line), locx is converted to a guide correction (GCx) using: 



GCx=<^ 



{locx - 40) when locx > 20 

locx when locx < 20 (8) 

when locx = 20. 



The simulations presented in Fig[T]were run 10000 times, and all apart from 7 (0.07%) predicted 
the correct guide correction required to align the reference and comparison projection. Each of the 
7 erroneous predictions were incorrect by a single pixel only. To obtain sub-pixel resolution for 
real CCD data, quadratic interpolation is performed across the three array elements surrounding 
^Xmax- This process is simultaneously carried out for the Y direction. We note that in the case 
of an extremely wide or distorted field a subsection of the CCD array centred on the target of 
interest would be be used to calculate guide corrections. Working in Fourier space allows quick 
cross correlation of the two projections created for each CCD axis. The calculation of the image 
shift is typically returned in a fraction of a second, see Tables [1] and [2] for a sample of off-sky guide 
corrections measured by DONUTS and their corresponding computation time. 



3. OfT-sky tests 

To demonstrate the robustness of DONUTS we tested it off-sky using archive data from four 
instruments, each with a different PSF, ranging from heavily defocused (FWHM ~ 40 pixels) 
to under sampled (FWHM < 2 pixels). The entire CCD window was used to calculate guide 
corrections, therefore using the positional information from all stars in the field, rather than from 
one star only. The remainder of this section describes briefly the instruments, observations and 
results from our off-sky testing of DONUTS. Since raw images are used for determining guide 
corrections in on-sky implementations, we have not reduced the archival images used here. To 
investigate the effects of reducing images before analysis with DONUTS, we bias subtracted, dark 
current corrected and flat fielded the test data described in Section 13.41 We find no improvement 
in the residual errors from DONUTS using reduced data. Assuming we have previously created 
a master bias, dark and flat field image, the reduction of each science image before analysis with 



DONUTS increases the calculation times in Tables [T] and [2] by at least a factor of 100. Therefore 
our assumption of analysing raw data only is justified. 

As in initial test, a random image from each instrument described below was chosen as a 
reference frame and subsequently shifted manually in the X and Y directions 200 times using IRAFg. 
The shift magnitudes were chosen randomly between ±40 pixels in each direction. DONUTS was 
then run on the series of manually shifted images to recover the known shifts applied. A sample of 
the results from this test for each instrument is given in Table [TJ See Figs [2][5] for examples of the 
PSFs tested off-sky. 

As a second, more realistic, test DONUTS was run once on the complete data series of each 
instrument below, measuring the required guide correction with respect to the first image in the 
set. This test is more realistic as the overall intensity of each individual image varies throughout 
each data set due to changes in airmass and seeing. The images were then shifted by the required 
correction as determined by DONUTS using IRAF. DONUTS was then run a second time on 
the corrected data to quantify the residual error. A sample of the guide corrections and their 
corresponding residuals is given in Table [2j The effects of changes in seeing on DONUTS is 
investigated in Section | 



3.1. Isaac Newton Telescope - Wide Field Camera 



A transit of the Super WASP dPohacco et ahlboOfil ;) exoplanet candidate ISWASP J161732.90+242119.0 



was observed on 2010-05-30 on the Isaac Newton Telescope (INT), with the Wide Field Camera 
(WFC). To increase the signal-to- noise ratio per unit time and to alleviate saturation the telescope 
was heavily defocused to ^ 40 pixels (13'.'2) FWHM (see Fig [2]). To increase the cadence, the WFC 
was also windowed from 4k x 2k to 2285 x 1465 pixels. Defocusing the WFC also defocuses the 
autoguider as the CCDs are in the same focal plane, and so the observations were guided manually, 
re-centering the telescope pointing every few minutes. A total of 167 manually guided images were 
obtained and subjected to the two tests described above. Regions of the CCD containing dead 
columns were excluded from guiding calculations. Small samples of the off-sky guide measurements 
and their residuals are given in Tables [T] and El 

The average residual errors from DONUTS when shifting a random image manually 200 times 
were only 0.01 pixel in each of the X and Y directions. The average residuals measured from the 
complete data series of 167 images were 0.02 and 0.01 pixels in the X and Y directions, respectively. 
The pixel size of the WFC is 0'.'33 pixel"^, therefore the theoretical average error from DONUTS 
corresponds to 0'.'007 and 0'.'003 on-sky. This, of course, is much smaller than the finest control of 
the INT (~ O'.'l) and the real error from DONUTS will most likely be dominated by the telescope's 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



limitations to make small guide corrections. 



3.2. Next Generation Transit Survey - Prototype Telescope 

On 2010-06-29 a total of 1 50 images with exposure times of 60 s were obtained of the transit 



ing super-Earth host GJ 1214 (jCharbonneau et alJl2009l ) during photometric testing of the NGTS 



prototype system (hereafter NGTS-P) in 2009/10. NGTS-P suffered from mechanical flexure be- 
tween the separate science and autoguiding telescopes as well as imperfect polar alignment. This 
was seen as a slow drift in stellar positions across the COD in the Y direction with a slight drift 
in X (3.79 pixel h"^ total). The telescope was defocused to 4.8 pixels (25'.'7) FWHM (see Fig [3]) 
to minimise the effects of the drift. The observations of GJ 1214 were analysed with DONUTS 
using the two methods previously described; small subsets of the off-sky guide corrections and their 
residuals are given in Tables [1] and [21 

The average residual errors from DONUTS when shifting a random image manually 200 times 
were 0.16 and 0.11 pixels in the X and Y directions, respectively. The average residuals over 
the complete data set of 150 images were 0.12 and 0.08 pixels, also in the X and Y directions, 
respectively. The plate scale of NGTS-P is 5'.'36 pixel"^ resulting in on-sky residuals of 0'.'643 and 
0'.'429 in X and Y, respectively. As in Section 13. H the average residuals from DONUTS off-sky 
testing at NGTS-P are also close to the limit of the telescope's repointing capability. Therefore, 
as before we suspect that the real error in autoguiding with DONUTS will be dominated by the 
response of the telescope. 



3.3. William Herschel Telescope - AG2 

AG2 was a temporary instrument with a frame-transfer CCD mounted at the folded Cassegrain 
focus of the Wilham Herschel Telescope (WHT). It had a Field of View (FOV) of 3f3 x 3f3 and a 
plate scale of 0^.^4 pixe l" ^ and was used for transit timing variation surveys of known exoplanets (e.g. 



Hrudkovaet al.ll20ld ). On the night of 2007-09-19 a transit of WASP lb (|Collier Cameron et al 
20071 ) was observed with 10 s exposures. The field of AG2 was significantly distorted resulting in an 



asymmetric PSF with a mean FWHM of 5.71 pixels (see FigH]). AG2 images were cropped before 
off-sky testing as they contained many pixels that were located outside the beam of the telescope 
and hence not illuminated. Observations of WASP lb were classically autoguided with the WHT's 
Gassegrain autoguider and hence showed very little movement in the stellar positions. Also, due to 
the small FOV of AG2, only 3 objects were visible on the CCD. The data were analysed using the 
methods previously described and samples of the results are given in Tables [1] and El 

The average residual errors from DONUTS when shifting a random image manually 200 times 
were 0.006 and 0.005 pixels in the X and Y directions, respectively. The average residuals from the 
raw data set were 0.002 and 0.006 pixels in the X and Y directions, respectively. This was expected 



as the stellar positions on the CCD were quite stable over time. As before, DONUTS appears to 
perform remarkably well regardless of the instrumental PSF, predicting guide corrections accurately 
at the sub-pixel level. 



3.4. Next Generation Transit Survey 

In August 2012 the first telescope (hereafter NGTS-1) for the full NGTS facility was delivered 
to the Geneva Observatory, Switzerland for characterisation and testing before being shipped to 
Paranal, Chile. The optical system of NGTS employs 20 cm aperture telescopes with fast focal 
ratios of f/2.8 made by Astro Syteme Austria (ASA) and Andor Technologies iKon-L, 2k x 2k 
CCD cameras with deep-depletion. Each telescope has a FOV and plate scale of 2? 83 x 2? 83 and 
4'.'96 pixel"^, respectively. On 2012-08-17 a series of 600 images with 4 s exposure times was taken 
during tests of the telescope's tracking performance. The images were focused such that the stars 
became under sampled with FWHM < 2 pixels (see Fig [5]). It was found that after precise polar 
alignment the system was remarkably stable and the stars showed very little long term drifting 
(< 1 pixel h~^) when tracking only. Analysis of NGTS-1 data with DONUTS off-sky using the 
methods described above also showed that the system required little correction, see Table [2] for a 
sample of the results. 

The average residual errors from DONUTS when shifting a random image manually 200 times 
were 0.18 and 0.11 pixels in the X and Y directions, respectively. The average guide corrections 
required to recenter the raw data set were 0.29 and 0.36 pixels in the X and Y directions, respectively 
with average residuals of 0.09 and 0.08 pixels in X and Y over the 600 corrected images analysed. 
The average residuals correspond to on-sky movements of ~ 0'.'40 which is slightly larger than for 
the instruments previously described due to the larger pixel size, but still close to the limit of the 
mount's pointing ability. Several individual residual measurements (Xr) on the x-axis of NGTS-1 
can be seen in Table [D to be slightly elevated. This was caused by an incorrect setup of the CCD 
clocking voltages which induced streaking of saturated stars in many of the images analysed. This 
issue has since been resolved and optimisation of DONUTS for the full NGTS facility is ongoing. 



3.5. Limiting signal-to-noise ratio and V magnitude of a single guide star 

When autoguiding with DONUTS it is important to quantify the limiting signal-to-noise ratio 
of the guide stars (including CCD noise) below which guide corrections become unreliable and 
therefore the algorithm breaks down. To ascertain this worst case scenario we created a series 
of 50 X 50 pixels image stamps, each containing a single star in focus, with a Gaussian profile 
(FWHM = 2.5 pixels). The signal-to-noise ratio (SNR) of the star was varied between 0.5 < SNR < 
20.0 in steps of 0.5. At each SNR the image of the artificial star was shifted manually 200 times in 
the X and Y directions. The shifts in each direction were randomly selected from values ranging 



between —10 and 10 pixels, with 0.01 pixel resolution. The same set of 200 shifts were used for all 
SNRs tested. 

We define a stringent, maximum residual error and standard deviation in the residual shifts 
from DONUTS of 0.3 pixel. This level of precision was surpassed when the SNR ratio of the 
manually shifted artificial star reached SNR ~ 9. Therefore, we place a conservative limiting 
SNR = 15 for the operation DONUTS. 

In Fig [6] we investigate as a function of exposure time the limiting V magnitude of a single 
guide star needed to give SNR = 15. We calculate contours for a 1 m telescope aperture (these 
can be scaled to different apertures) for bright and dark sky conditions, and for a range of CCD 
parameters (see Table [3]). We assume a baseline model CCD with a read noise of 10 e~ and dark 
current of 0.2 e~ s~^ pixel"^ and investigate the effects of increasing these noise parameters to 
simulate lower quality CCD cameras. Using Fig [6] and, 

we can calculate a range of allowable exposure times (Texp) and guide star brightnesses, where Tc 
is the maximum allowable dwell time before the pointing of the telescope needs corrected and Tr is 
the read out time of the CCD. Assuming Tc = 120 s and Tf = 20 s the exposure time must satisfy 
^exp ^ 100 s, resulting in limiting guide star brightnesses of y ~ 19 and 20.5 for our baseline model 
during bright and dark sky conditions, respectively. 

CCDs with higher read noise require brighter guide stars to reach the required SNR=15 when 
observing at high cadence, this is because the read noise from the CCD is significant in comparison 
to the noise from the guide star and sky background. As the exposure time is increased, the dashed 
lines of increased read noise in Fig [6] tend to converge with the baseline model, this convergence is 
most pronounced during full moon as the noise from the sky quickly begins to dominate. CCDs with 
increased dark current require slightly brighter guide stars when observing with long exposures as 
the flux from dark current builds up quickly. We can see this as a very slight divergence of the dotted 
lines in Fig [6] from the baseline model. However, by far the most dominant factor in determining 
the limiting magnitude of a guide star, for a given exposure time, is the sky brightness. Figure 
[6] shows a change in limiting V magnitude of approximately 2 mag between full and new moon. 
Defocused stars were not analysed in these simulations as by definition defocusing is only applied 
when the target is bright and therefore does not represent the limiting case of autoguiding low SNR 
observations with DONUTS. 



3.6. Changes in seeing 

When observing a field of heavily defocused stars, changes in seeing have very little effect on 
the effective PSF crpsF as the focus term in 
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0-pSF = y f^seeing + C^focus' (10) 

tends to dominate, where cJseeing and cifocus are the uncorrelated FWHMs of the seeing disk and 
telescope focus in arcsec, respectively. However, when observing close to focus, changes in seeing 
may dramatically alter the PSF. This is most apparent for instruments whose pixel scale is less 
than the typical site seeing for the observatory where it is located. 

To test the effects of changes in seeing on DONUTS we conducted the following simulation. 
We assume that a reference image for a given field has been observed in focus and under good seeing 
conditions, however on a subsequent observation of the same field the seeing has deteriorated. To 
simulate the worst case scenario our reference image was blurred with a Gaussian blurring filter 
and shifted manually 200 times. As before, the shifts were chosen from a range of ±40 pixels in 
both the X and Y directions. This process was then repeated for 8 levels of seeing deterioration, 
ranging from a 0% to 125% increase in the average stellar FWHM. 

We applied this method to our under sampled NGTS-1 data, using the same reference image 
as used in Section 13.41 The same 200 shifts applied during our off-sky testing were applied to the 
reference frame and all but the reference image were blurred using a Gaussian blurring filter to the 
FWHMs given in Table HI In a similar manner to the first off-sky test in Section [31 DONUTS was 
then run on the shifted and blurred data and the measured shifts were compared to the known 
offsets applied. The average residuals for each level of simulated seeing change are compared to the 
normal data in Table HI As can be seen in Table H] changes in seeing have only a marginal effect on 
the residual errors of DONUTS. 



4. On-sky test 



DONUTS was first tested on-sky at the NITES telescope (JMcCormac et al-llAfcA submitted ) 



on La Palma. The NITES telescope is a remotely operated 0.4 m Schmidt-Cassegrain telescope 
with a 1024 x 1024 pixel CCD camera with deep depletion technology. The FOV and plate scale 
of the telescope are 11'26 x 11'26 and 0'.'66 pixel"^. The NITES telescope suffers from periodic 
error (P ~ 4 min) in right ascension caused by mechanical inaccuracies in the drive system. During 
commissioning the periodic error was modelled and corrected, reducing it from approximately 
±30 to ±6 pixels. As DONUTS applies guide corrections between science images the periodic 
error essentially limits the maximum exposure time of NITES to tmax = 30 s without the use of 
conventional autoguiding. Figure [7] shows the effects of enabling DONUTS at the NITES telescope. 
Figure [7| (X symbols) shows the drift in position of a star observed on 2011-02-18 with DONUTS 
disabled. The observations lasted 3.5 h and a drift can clearly be seen in declination. This is caused 
by a slight polar misalignment of the telescope. The vertical spread in the drift profile is caused by 
the periodic error of the mount. Figure [7] (cross symbols, see subplot insert) shows data taken over 
4.3 h on 2012-11-19 with DONUTS enabled. A clear improvement over tracking only can be seen. 
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The higher the cadence of observations made with the NITES telescope the more times the periodic 
error is sampled and corrected by DONUTS. If deployed on a mount with little or no periodic error 
(e.g. NGTS-1) we expect extremely stable tracking performance over long periods of time. 



4.1. Discussion and summary 

Our initial simulations of DONUTS in Section [2l using a simulated star which had been slightly 
defocused returned encouraging results. This test showed that DONUTS is capable of measuring 
translational shifts between two 1-D projections accurately. Providing the limiting SNR = 15 is 
achieved, changing the overall intensity of the comparison projection in Fig [1] has no effect on 
determining the shift correctly. Also, subsequently tweaking the individual pixels of the projection 
randomly by a moderate amount, and varying the overall intensity of the comparison frame has no 
effect on the shift measurement. The simulations in Section [2] were run 10 000 times, with 0.07% 
of the measurements being incorrect by a single pixel only. 

In Section [3] we demonstrated DONUTS capability of rapidly calculating guide corrections 
using a full CCD array and a variety of PSFs, while maintaining an average residual error well 
within the sub-pixel (< 0.2 pixel) regime. We define a conservative limiting SNR in the worst 
case scenario of SNR = 15 and investigate the effects of changing sky brightness, dark current 
and read noise on the limiting V magnitude of a single guide star in the general case of a 1 
m telescope. The benefits of spatial stability combined with the ability to eliminate mechanical 
fiexure, while simultaneously defocusing (if desired) will aid astronomers in pushing the limits of 
photometric accuracy to the sub-millimag level in many areas of astronomy. For example, near 
infrared observations of faint brown dwarfs or secondary eclipses of transiting exoplanets, which 
are typically carried out in focus to minimise the eff'ects of the infrared background, could benefit 
as much from DONUTS as heavily defocused optical observations of bright transiting exoplanet 
host stars. 

As seen in Table [2] several guide corrections with magnitudes < 0.1 pixels become difficult to 
correct, however at this level of spatial stability we are theoretically already able to meet our goals 
for DONUTS in providing sub-pixel stability (< 0.2 pixels) for high-cadence time-series photometry. 

Section d] clearly highlights the main limitation of the DONUTS autoguiding algorithm. De- 
ployment of DONUTS on a poorer quality mount results in a general localisation of stellar positions 
rather than the absolute sub-pixel fixing we desire. However, as shown in Section [3.41 the NGTS- 
1 telescope is extremely stable and we therefore expect unprecedented tracking performance for 
NGTS when DONUTS is enabled. 

We note that, over wide FOVs, spatial differential refraction caused by the atmosphere may 
introduce a slight error if the entire FOV is used to calculate guide corrections. Under such 
conditions the guide corrections may be calculated using a small subsection of the CCD surrounding 
the target of interest. Our investigation into the effects of changes in seeing on DONUTS also 
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returned promising results. The average PSF in the comparison frames tested was increased by up 
to 125% with almost no effect on the level of residual error from DONUTS. This was expected as 
the maximum correlation between two Gaussian profiles occurs when their peaks overlap, regardless 
of the increase in FWHM. As so long as the telescope is not physically defocused by a large amount 
(such that the secondary obstruction is visible) between acquiring the reference and comparison 
images we expect DONUTS to handle slight variations in PSF and overall image intensity in a 
robust manner. 

We have presented a new autoguiding algorithm designed to correct a telescope's pointing 
between the science images of high-cadence time-series photometry, with the goal of producing 
stable telescope tracking at the sub-pixel (< 0.2 pixel) level. Our tests have shown that autoguiding 
with DONUTS will theoretically be superior than conventional autoguiding up to an exposure 
time where telescope tracking errors become a significant fraction of the PSF size. DONUTS 
simultaneously eliminates mechanical flexure issues and atmospheric refraction effects as a target 
is tracked through a range of airmasses. We conclude that the algorithm is significantly robust 
at handling different PSF sizes and shapes and is robust against changes in seeing. We find an 
average residual error < 0.2 pixels over the range of PSFs tested. The main limitation of DONUTS 
is expected to be the response of a given telescope to small guiding corrections. One possible way 
of circumventing this issue is to characterise the telescope response and implement a control loop, 
e.g. PID, to account for over or under correction made by the telescope. Tailored versions of 
DONUTS are currently under development for the NGTS facility at Paranal in Chile and for the 
1 m SuperWASP Qatar Telescope (SQT) at the Observatorio del Roque de los Muchachos, on the 
island of La Palma in the Canary Islands. 

Facilities: NGTS, NGTS-P, ING:Newton (WFC), ING:Herschel (AG2). 
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Table 1. A small sample of the shifts applied manually to a reference image from each of the 4 instruments tested. Xm and 
Yjn are the manually applied shifts in pixels. Xj-m and Yj-m are the differences in pixels between the shift applied and that 

measured by DONUTS and teak is the time in ms of the calculation. 
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Table 2. A small sample of predicted guide corrections calculated off-sky for the 4 instruments tested. X and Y are the 

required guide corrections in pixels between the reference image and a random sample of 10 images from the observing block. 

Xj. and Yr are the residuals in pixels measured by DONUTS after the corrections have been applied manually in IRAF and 

tcaic is the time in ms of the calculation. 
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Table 3. Assumed parameters used for calculating the limiting V magnitude in the worst case 
scenario of a single, low signal-to-noise guide star. 



Parameter 



Assumed Value 



Limiting SNR 

Telescope aperture 

Camera platescale 

Read noise 

Dark current 

FWHM star 

Photometry aperture radius 

Total efficiency'' 

Vega flux 

Sky level at full moon'^ 
Sky level at new moon'' 



15 

1.0 m 

0.5 arcsec pixel"'^ 

10 e' 

0.2 e~ s^i pixcl^l 

2.5 pixels 

2.5 pixels 
50 % 

9.6 X 10^'' photons s"! 
0.555 iim 

0.089 lira 

18 mag arcsec"^ 

22 mag arcsec"^ 



fira 



''Includes atmospheric, telescope and filter throughputs plus CCD 
quantum efficiency in the V band. 

''We assume the sky brightness levels for La Palma listed at 
|http://www. ing.iac.es/PR/newsletter/news6/tell. html 
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Table 4. Simulated seeing changes after applying a Gaussian blurring filter. Row 1 represents 

the original data, X^ and Yj- are the average residuals for each 200 image series at each simulated 

seeing level. 5fwhm is the percentage increase in FWHM at each stage with respect to the 

original data (row 1). 



FWHM Xr Yr 5pwHM 

(pixels) (pixels) (pixels) % 



1.32 


0.18 


0.11 





1.58 


0.18 


0.10 


20 


1.78 


0.18 


0.10 


35 


1.98 


0.19 


0.10 


50 


2.18 


0.20 


0.09 


65 


2.37 


0.20 


0.09 


80 


2.58 


0.21 


0.09 


95 


2.78 


0.22 


0.09 


110 


2.98 


0.23 


0.09 


125 
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Fig. 1. — Simulated guide star testing with a slight defocus, such that the simulated PSF becomes 
double peaked. The solid red line represents a 1-D projection of a typical guide star in the X 
direction, where a slight telescope defocus has been applied. The solid blue line represents the same 
guide star with a uniform change in intensity applied to all pixel elements. The solid green line 
represents the same as the solid blue line plus a random pixel-by-pixel change in intensity where 
the random values are selected using a Gaussian distribution with a mean of and a standard 
deviation of 1 from ±20% of the intensity modified (blue) pixel value. The dashed blue and green 
lines represent the IFFT, (px, from Eq[7]of the solid blue and green comparison projections above. 
The location, locx, of the peak c/iXmax! ™ 'Px, with respect to the x-axis origin is converted to a guide 
correction using Eq [HI The same treatment is also applied to the Y axis projections to calculate 
the guide correction in the Y direction. 
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Fig. 2. — Example PSF from the INT/WFC used in off-sky testing. 
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Fig. 3. — Example PSF from the NGTS-P telescope used in off-sky testing. 
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Fig. 4. — Example PSF from the WHT/AG2 used in off-sky testing. 
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Fig. 5. — Example PSF from the NGTS-1 telescope used in off-sky testing. 
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Fig. 6. — Limiting V magnitude of the DONUTS autoguiding algorithm as a function of exposure 
time for the worst case scenario where only one guide star is available. The solid lines show our 
baseline model of a 1 m telescope and assuming the parameters given in Table [3l The dashed lines 
represent an increase in read noise of 10 e~ pixel^^ on top of the baseline model (20 e~ pixel"^ 
total) and the dotted lines simulate a considerable dark current increase from 0.2 to 1 e~ s~^ 
pixel"^. The red and black contours represent full and new moon observations, respectively and 
the insert shows a zoom around short exposure times. As expected, increasing the read noise has 
a strong effect on the limiting guide star brightness at short exposure times. Increasing the level 
of dark current has a noticeable effect on the limiting guide star brightness at long exposure times, 
especially during new moon conditions. Increasing the CCD noise during full moon conditions has 
little effect when Texp > 30 s as the noise from the sky begins to dominate. 
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Fig. 7. — Example plots of stellar positions across the NITES CCD with DONUTS disabled (X 
symbols) and enabled (cross symbols, see subplot insert). The position of each star has been 
normalised to (0,0) using its position in the first image. Although the NITES telescope suffers 
badly from periodic error in the RA axis there is a clear improvement in using DONUTS over 
tracking only. 



